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Abstract 

(N 

We compare different treatments of the constraints in canonical 

■ quantum gravity. The standard approach on the superspace of 3- 
. geometries treats the constraints as the sole carriers of the dynamic 

content of the theory, thus rendering the traditional dynamical equa- 
£f~^ ' tions obsolete. Quantization of the constraints in both the Dirac and 

■ ADM square root Hamiltonian approaches leads to the well known 
Q ! problems of time evolution. These problems of time are of both an in- 

terpretational and technical nature. In contrast, the geometrodynamic 
^jy quantization procedure on the superspace of the true dynamical vari- 

ables separates the issues of quantization from the enforcement of the 
constraints. The resulting theory takes into account states that are 
^ ■ off-shell with respect to the constraints, and thus avoids the problems 

^ . of time. We develop, for the first time, the geometrodynamic quan- 

tization formalism in a general setting and show that it retains all 
essential features previously illustrated in the context of homogeneous 
cosmologies. 
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1 Introduction. 



The standard approach to the canonical quantization of gravity P3I2| is based 
on the classical dynamic picture of the evolving 3-geometry of a slicing of a 
spacetime manifold. The slicing is essentially a reference foliation of the 
spacetime manifold (endowed with a 4-geometry) with respect to which 
canonical variables are assigned. It is usually parametrized by a time co- 
ordinate t and tied to the enveloping spacetime by the lapse function N and 
the shift vector iV\ The canonical variables are the 3-metric gik on a spatial 
slices of the foliation induced by the spacetime 4-metric, and the matrix -K lk 
of their canonically conjugate momenta. The latter is related to the extrinsic 
curvature of S when it is considered to be embedded in the spacetime. 

The customary variational procedure, applied to the Hilbert action ex- 
pressed in terms of the canonical variables, produces Hamilton equations 
which describe the time evolution of the canonical variables. The Hamilto- 
nian is given by NTi + N{H % , where 7i and H? are functions of the canon- 
ical variables and their spatial derivatives. The procedure is not extended 
to the derivation of the Hamilton-Jacobi equation in the usual manner, as 
such an equation with the chosen set of canonical variables appears to be- 
come meaningless when the general covariance of the theory is implemented 
(cf. section for a more precise description). 

Alternatively, general covariance can be introduced in the variational 
principle from the outset, by requiring that the action is invariant with 
respect to variations of the lapse and shift. This leads to the constraint 
equations (to simplify notations, we omit indices on components of g and tt 
in all equations in this section) 



which are imposed on the canonical variables of each slice. An important 
feature of general relativity is that its dynamics is fully constrained. It can 
be shown that if the geometry of spacetime is such that the constraints are 
satisfied on all slices of all spatial foliations of spacetime, then the canonical 
variables necessarily satisfy the Hamilton evolution equations. This feature 
is often referred to as a key property of general relativityjT] and is used to 
argue that the entire theory is encoded in the constraints, with the conclusion 
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that the Hamilton equations are redundant and can be ignored in dynamical 
considerations. Substitution of SS/Sg in the place of p in the constraint 
equations leads to a new set of equations 
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the first of which is considered to be the Hamilton-Jacobi equation (see sec- 
tion EJ). Arguments appealing to the variational principle on the superspace 
of 3-geometries support this assertion. Detailed arguments and the interpre- 
tation of other equations may be found elsewhere jl]. 

Dirac's procedure of canonical gravity quantization is based directly on 
this Hamilton-Jacobi equation, and produces a quantum theory that con- 
sists of the Wheeler-DeWitt equation together with commutation relations 
imposed on all canonical variables. 

The ADM square root quantization procedure is also based entirely on 
constraints, but in this procedure the set of canonical variables is split into 
two subsets: the embedding variables (four of them altogether; one slicing 
parameter Q and three coordinatization parameters a) and the true dynam- 
ical variables (3 (two of them) jSJ U\ ■ The constraints are then solved with 
respect to the momenta conjugate to the embedding variables. After substi- 
tuting 5S/5Q, 5S/5a for p n , p a , (where S is the principal Hamilton function) 
one of the resulting equations (the equation for the momentum conjugate to 
the slicing parameter) is identified with the Hamilton-Jacobi equation, and 
its right hand side yields an expression for a new (square root) Hamiltonian. 
The quantization is based on this equation, yielding a quantum theory that 
consists of the Schrodinger equation and commutation relations imposed on 
the true dynamical variables and their conjugate momenta. 

Describing the time evolution of quantized gravitational fields is extremely 
troublesome in both formulations. Any attempt to introduce time that can 
be used in a way similar to time in quantum mechanics, or in quantum field 
theory on a flat background, invariably leads to the notorious problems of 
timejS]. Attempts to introduce time in a universal way externally (such as 
the readings of specially designed clock) have been unsuccessful. There are 
reasons to believe that this is impossible pQ, whether the clock is believed 
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to be gravitationally defined (i.e. the readings depend only on the variables 
describing the gravitational field), or a matter clock (the readings depend on 
both gravitational and matter variables). 

The conceptual difficulties (such as the problem of functional evolution 
and the multiple choice problem, in Kuchaf's terminology) emerge due to 
the dual nature of time parametrization in general relativity. If spacetime is 
considered as a manifold, it can be coordinatized and sliced in an arbitrary 
manner. However, this is not sufficient for the description of geometrody- 
namic evolution, since both slicing and coordinatization need to be tied to the 
metric of spacetime. The standard approach in classical geometrodynamics 
involves lapse and shift. These quantites determine the slicing of spacetime 
(i.e. the slicing condition) and its coordinatization. On a physical basis the 
lapse (the slicing) represents the reading of the clock of the observer whose 
worldline is perpindicular to the spacelike hypersurface. While the shift rep- 
resents the displacement of the spatial coordinates in time away from such an 
observer. Here, the dynamics is determined solely interms of the constraints 
only if there exists a unique spacetime metric. This metric might not be 
known until the geometrodynamic problem is solved. While its existence is 
not a problem in classical general relativity, there is, in general, no possibility 
of assigning such a unique metric of spacetime in canonical quantum gravity. 

Given these problems, continuing to use the constraint equations as the 
foundation of geometrodynamics becomes problematic, to say the least. Quan- 
tization of the dynamical picture based on the constraints is essentially equiv- 
alent to restricting the states of the resulting quantum systems to a "shell" , 
which is determined by constraints that are classical in origin. An attempt 
to undertake a similar action in quantum mechanics or quantum field theory 
would be quite disastrous under all but very carefully selected conditions. 

One way to resolve this dilemma would be to weaken the requirement 
of covariance, essentially discarding it in dynamical considerations and re- 
covering it by imposing symmetries on solutions only to the extent and in 
the sense that is allowed by dynamics. General covariance in the traditional 
sense should be recovered in the classical limit. This requirement should 
determine, at least partially, what constitutes the classical limit of quantum 
geometrodynamics. 

This can be achieved if York's analysis of gravitational degrees of freedom[H] 
is taken into account and actively implemented. The Hamilton-Jacobi equa- 
tion takes its traditional form, familiar from classical mechanics or electrody- 
namics. The resulting description remains equivalent to the one commonly 
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accepted in classical geometrodynamics. However, quantization based on 
this new Hamilton-Jacobi equation provides an appropriate interpretation 
of the conceptual problems of time, making them quite natural statements 
concerning the properties of gravity quantization. It also seems to avoid the 
technical problems of time, such as the Hilbert space and spectral analysis 
problems, as it produces a Schrodinger equation for the state evolution, and 
the Hamiltonian does not include the square root operation. The proce- 
dure has been described elsewhere in the context of homogeneous quantum 
cosmologies [HI Ej- 

In this setting time can be introduced as a slicing parametrization on a 
spacetime manifold, and tied to the metric structure, without contradiction. 
The metric interpretation of time is coupled with geometrodynamic evolu- 
tion, while the true meaning of time becomes completely determined only 
after the geometrodynamic evolution problem has been solved. In a sense, 
the quantum geometrodynamic configuration and time emerge together, and 
the meaning of the clock readings is influenced by the quantum gravitational 
system j7]. 

As noted above, the geometrodynamic approach differs from the standard 
one in its treatment of the constraints. This difference is rather subtle on the 
classical level and does not result in different predictions. On the quantum 
level, however, it becomes fundamental. 

In previous work we have considered different aspects of geometrodynamic 
quantization in particular cases of homogeneous quantum cosmologies. Here, 
we develop, for the first time, the geometrodynamic quantization formalism 
in a general setting and show that it retains all of the essential features 
previously illustrated in the context of homogeneous cosmologies. 

To achieve transparency of the discussion and to provide an appropriate 
platform for future applications, we start by reviewing the issue in classical 
geometrodynamics. We begin from the Lagrangian formulation, and then 
describe the transition to the Hamilton and Hamilton-Jacobi equations. 

2 Constraints on the Configuration Space of 
3-Metrics. 

We start from the 3+1 Lagrangian expression for the action, obtained from 
the standard Hilbert action by expressing the 4-metric in terms of lapse N, 
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shift Ni, and 3-metric g^. After eliminating total time derivatives and total 
divergences jl], the action becomes 



16n 



R — (TrK) 2 + Tr(K 2 ) N^/g~d 3 xdt= Cd 3 xdt 



where 



1 



16n 



R — (TrK) 2 + Tr(K 2 ) Ny/g~. 



(5) 



(6) 



Following standard convention^) we drop the factor of 1/167T from the grav- 
itational action in the remainder of this paper. In vacuum the equations are 
equivalent, and the factor can be trivially included when matter is present. 
In the equations above K is used as a shorthand notation for the extrinsic 



curvature Kij, defined as 



1 

2iV 



N i{j + N. 



9ij 



(7) 



R is the 3-curvature associated with the 3-metric gij, g is the determinant 
of this 3-metric, and g^ = 

Components of the 3-metric g^ are treated as dynamical variables or 
(functional) coordinates of the geometrodynamic configuration space. Fol- 
lowing the standard prescription for transitioning from the Lagrangian de- 
scription of dynamics to the Hamiltonian one, we introduce momenta 7r u 
conjugate to the dynamical variables g^, with 
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dC 

dgij 



(8) 



Computing the right hand side of this expression involves determining the 
derivatives of K = Tr(K) = g^K^ and of Tr(K 2 ) = g im gi k K mj K kl with 
respect to the g n \ using 



dK, 



dg k r 



— s k s m - 

2N 1 r 



This computation yields 
<9TrK 



dg 



nl 



2N 



2N 



(9) 



(10) 
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<9(TrK 2 
dg, 



1 



id 



2N 
1 



# g J d m d j K ki + g g J K mj d k o j 



2N 



in Ik 



K kl + g^g lm K n 



-— K 



N 



g g ^ki + g J g - 

which results in the expression for the momenta ir^ defined by 

* - y/g (Kg* - 



7T 



(11) 



(12) 



In what follows we use notations II = 7r* J and n = Trll. The last equation 
implies 



n 



^g(3K-K) = 2 v fgK 



K 11 = -Fig 13 

V9 V2 



7T 



*3 



TYK 



if™ 2 - — 
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(13) 
(14) 

(15) 



which allows one to express £ as a function of gij and tc* only, giving 

i2 



C = N 



g%R-g-* (H_-Trn 2 



The standard definition of the Hamiltonian 

Hwdw = n v 9ij — C (Mod total divergence) 



(16) 



(17) 



can be transformed to express H-wdw in terms of g^ and the conjugate 
momenta. In order to achieve this, we use 



gij = % + Nj\i - 2NK ij 
which allows us to write the first term of (JT7j) as 

"'•'.'/,; = (% + - 2N^K ij = 2N^ - 2N^K ij 



2 (N t n% - 2Ni (tt% - 2.X~"K, j . 



(18) 



(19) 



The first term in this expression is the total covariant divergence of a vector 
density It can be expressed as the total divergence of a vector field, 



(iViTT*),. = ^g [Ni{KgV - K*)]^ = [y/gNiiKg** - K ij 



(20) 



and removed from the final expression for the Hamiltonian. The third term, 
after substituting as given by (|14|). reduces to 

2Nn ij K ij = 2g~*N - TrU^j , (21) 

and the expression for TLwdw takes the form 



n WDW = N l (-2^) +N 



g~* hiU 2 -^-) -g*R 



(22) 



It is common practice to write the Hamiltonian Hwdw hi the form 

H WDW = NH + N l H l (23) 

where 

H = a -\{ TrTT 2 - ■ 2 
is called the superhamiltonian, and 



H = g-^ (ttU 2 -^-) -g$R (24) 



W = -2(t? 3 ) (25) 

are the supermomenta. The action in Hamiltonian form can thus be written 

as 

I c = J (jr^gij - Hwdw) d 3 x dt 

= J (tt%- -NH-NiTi?) d 3 x dt. (26) 

Variations of this action with respect to tt 1 ^ and gij produces the Hamilton 
equations 

9Hwdw 

= -!™*™L, (28) 

while variations of the shift and lapse yield the superhamiltonian and super- 
momenta constraint equations 

W(<fc,7T<0 = (29) 
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W*(0y,7r«) = O, (30) 

which will be discussed in more detail after we develop convenient notations 
and formal machinery to handle the relevant questions. It is important that 
these equations (|27|) can be obtained by inverting the kinematic relations 
(|12|): therefore, they are independent of the Hamiltonian dynamics. 



3 Transformation of Variables. 

It is often convenient to replace with another set of variables (functions) 
which form the configuration space. For instance, in general analysis of the 
initial data problem or gravitational degrees of freedom the variables are split 
into the true dynamical variables, the scale factor and the gravitomagnetic 
vector (three components) that cannot be identified with components of the 
3-metric gij. A similar parametrization is used in setting up the problems of 
homogeneous cosmologies jU]. In general, the are expressed as 

9ij = 9ij(QA) (31) 

or 

9%j = gijiUiX 1 ) (32) 

where q A = Qa( x \ t) are assumed to be independent, and A = 1, . . . ,n q with 
n q < 6. In the generic case n q = 6 and all components of shift are present in 
the spacetime metric. If this is not the case, some symmetries have been used 
to fix the form of the metric, which typically leads to the loss of covariance. 
Some of the supermomenta constraints might be lost or not be independent. 
It is hard to list and consider all possibilities, but as a rule these degenerate 
cases do not cause essential troubles in any particular case. In what follows 
we will be interested mostly in the generic, non-degenerate case, although 
most of the conclusions will also be true for degenerate cases. 

The transition from the Lagrangian to the Hamiltonian action begins with 
the Lagrangian (JHJ). The extrinsic curvature is given by (J7J), except now g^ 
is given by 

g ij = ^ L q A = M A ij q A (33) 
dq A 

which results in the expression for K. 



Kij ~ 2N 



13 

1 



Niu + Nju - M A ijqA (34) 
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and 



dK 



1.7 



dq A 



— M A ih 
2N lv 



(35) 



where M A ij is a 3 x 3 symmetric matrix for each value of A. Alternatively, 
one can narrow it down to six components of gij such that i < j (right 
upper triangle) and consider pairs (ij) as collective indices, in which case this 
matrix becomes the Jacobian of transformation between q A and gij,i < j. 
The assumption that the variables q A are independent implies that the rank 
of this matrix is equal to the number of q A (less or equal to six). In the generic 
cases this rank is equal to six, which means that the system of equations (|33|) 
can be solved, providing expressions for q A in terms of g^. 
Equations (f33J) allow us to calculate the derivatives of 



TrK 

Tr(K 2 ) 

with respect to q A , which yields 



g im gi k K mj K kl 



<9TrK 

dq A 
g(Tr(K 2 )) 
dq A 



2N 



9 ij M A t3 



(36) 
(37) 



(38) 



M A mj 

2N 1 



K kl + g im g lk K,, 



rnj 



1 

"2N 



M A 



ki 



2N 

-—K ij M A 
N 



K^ m M A m j + K ik M A ki 



(39) 



The momentum p A conjugate to the variable q A is given by 



,.4 



P 



dC 

dq A 



^g- Kg'i.U ';,- K'i.U- 



M J 



(40) 



and comparison of this expression with (|T^|) yields the very useful relations 

P A = ^M A %r (41) 
, with K given by ()34j) . becomes a function of the new 



The Lagrangian £ 
variables 



C = £{q A , q A , N, Ni[ 



(42) 
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and the transition to the Hamiltonian action can proceed in the standard 
way, expressing the Hamiltonian 

H = p A q A - C, (43) 

in terms of canonical variables q A , p ■ This procedure presumably results in 
an expression for the Hamiltonian of the same type as above: 

U = NH(q A ,p A ) + NiH'iq^P^. (44) 

This is obvious in the generic case, when the rank of M A ij is maximal (six), 
since (I44J) can then be obtained by inverting ()31|) and (}4*Tj) and substituting 
the expressions for 7T 1 - 7 = -a %3 (q A , p A ) , together with (|3Tj). into PB^ - ipB]) . The 
resulting Hamiltonian 7i depends on the new variables q A ,p A , and inverting 
(jHIJ), followed by the substitution of q A = q A (gij) and (}4*Tj) into this new 
Hamiltonian, yields 

n(q A (g tJ ),p A (g tJ , tt"), N, N>) = H WDW ( 9ij , ^\ N, N t ) , (45) 

which implies 

n/r A ■ dHwDw dp A dH A &H , 

9 - = M - qA = ^^ = ^W A=M w (46) 

Thus 



and this system of equations has the unique solution 




- - w- (48) 

Just as in the previous section, these expressions for q A are equivalent to the 
definitions of momenta p A given by (|40|) . 

This logic fails when the number n q of variables q A is less than six. Never- 
theless, the basic structure remains the same. In particular, the Hamiltonian 
Tt (q A} p A , N, Nij can be expressed by ()44j) as a combination of the super- 
hamiltonian and supermomenta (neither of which depend on shift and lapse). 
Note, however, that this expression cannot be obtained by the simple inver- 
sions described above. Just as in the non-degenerate case, the Hamilton 
equations (J48|) remain equivalent to the definitions of momenta (J4(Jj) . 
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This procedure is more complex, and requires lengthier computations, 
but we believe that it leads to a greater understanding of the structure of 
the equations, and thus we present it here. As an additional benefit, this 
procedure provides explicit expressions for q^ and p A in terms of each other, 
as well as explicit expressions for the superhamiltonian and supermomenta. 

It is convenient to rewrite equation (|4(J|) for p A in the form 



V 



v'7/ (av-m/ \. ( - K L m A ij 



As before, we have 



K, 



i in 



2N 



Ni\ m + N m \i — M B i m q B 



which yields 

K — " lm 



9 lm K lm = — [g im (N l{m + N m \i) - g lm M B lm q B 



K ij = g il gl m K, 



tin 



2N 



9 a 9 ]m (Nn m + N m \i) - S V"* M B lm q B 



(49) 
(50) 

(51) 
(52) 



and substituting the last two expressions into (j4T?jl leads to 



V 



2N 



i9 U 9 jm -9 ij 9 lm ) M A l] (N l]m + N llm ) 
+{g il 9 jrn - 9 lJ 9 lm ) M\ M B lrn q B ] . (53) 
This can be written in a more compact form if we introduce the notations 

Qijlm = g il g jm _ g ij g lm ^ 

or its symmetrized version 



Qijlm _ _ g il g jm _|_ ^gjlgim _ gijg 1 ™ 



and 



Q AB = G ijlm M A ij M B lm 
With these notations, the momenta becomes 



p 



2N 



G ijlm M A.. + Ni]m) + qAB - b 

12 



(55) 
(56) 

(57) 



We have assumed that the q A are independent. This implies that the square 
matrix Q AB (with dimension n q x n q ) has rank n q and, consequently, is 
invertible. In what follows we use the same letters G and Q with lower 
indices for the elements of the inverse matrices, 

(G ijlm ) = (G^ lm y l (58) 
(Qab) = (Q AB Y\ (59) 

so that 

Gij kn G kn lm = 8\ 8 3 m (60) 

and 

Qab Q BC — 8 a- (61) 
These matrices have the quite obvious symmetry properties 

Qijlm Qjilm Qijml Qlmij (62) 

Gijlm Gjil m Gij m l Gl m ij (63) 

as well as 

Q AB = Q BA - Qab = Qba- (64) 

There are two more useful relations that follow from simple observations. By 
definition 

Qbc Q CA = Qbc M c i3 G ij lm M A lm = 5 A (65) 

which implies 

M B kn Q BC M c \ 3 & lm M A lm = 5 A M B kn = M A kn = 5 l k 8™ M A lm . (66) 

i &re a complete set of variables (g ik depend only on q A and do not 
have any other arguments) then this results in the two relations 

M B kn Q BC M c l3 G ij lm = 5{ 8™ (67) 

M B kn Qbc M c i3 - = G kni3 (68) 

that will be used, together with symmetries, in computations below. 

The expression (|57|) for the momenta can be considered as a system of 
linear equations for q B - Its solution 

q B = Qba G^ lm M A l3 (N l]m + N m]l ) + 2Ng^ Q BA p A (69) 
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expresses q B in terms of the momenta p A . 

We begin to calculate the Hamiltonian TC by expressing the extrinsic 
curvature Ki m , Tr(K) and Tr(K 2 ) in terms of momenta. The last term in 
square brackets in equation ()50|) for K\ m can be written in terms of momenta 
by substituting q B given by 



M B 



lm QB 



lm 



q BA G^ lm M\ (N k{n + N nlk ) + 2Ng~l Q BA p A M E \ 



lm 



fjk In 
I rn 



= Ni\ m + N m \i - 


Y2Ng-^Q BAV A M B lm . 


(70) 


This results in 






Kim = 


~9~^ Qba p A M B lm 


(71) 


TrK = 


-g- 1 iQ B AP A g lm M B lm 


(72) 


K ij = 


-g- 1 *QBAP A g u g jrn M B lm 


(73) 


Ka = 


-9^Qcdp d M c {j 


(74) 


which, together with Tr(K 2 ) 


= K l i Kij, allows us to express the Laj 


;;rangian 



C in terms of momenta: 
C = 



R — (TrK) 2 + Tr(K 2 )J y/g N 

g?R + g~$ Q BA Q CD G^ lm M c \ 3 M B 



lm P A P° 



= N 

= N [g^R + g^Q AB p A p B . 
In addition, (ffiHJ) implies 

p A q A = 2Ng~* Q AB p A p B + 2 Q AB p A M B lm G lm « 
which results in the expression for the Hamiltonian 
7i = p A q A — £ (Mod total divergence) 



(75) 
(76) 



N 



g- 1 *QABP A p B -g k *R]+2N ilj Q AB p A M B lm G lmi i. (77) 



The last term in this expression can be written as 

2 Ni\j Q AB p A M B i m G lrni i = (2N l Q AB p A M B lm G lmi i 



-N t {2Q AB p A M B lm G lm ^ 



(78) 
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The covariant divergence of a vector density is a total divergence and can be 
thrown out. This reduces the Hamiltonian to 



H = N 



g~* QabP A P B -g*R]+ N t (-2Q ABP A M B lm G lmi i) b (79) 
which is usually written as 

H = NHiq^p^+NtH'iq^pA) (80) 

where 

H(q A ,P A ) = <T* Qabp a p b -g*R (81) 

is the superhamiltonian, and 

n\q A ,p A ) = (—2 Qab p A M B i m Q lmi i\ (82) 



are supermomenta in the superspace of q A . It is trivial to verify, by computing 
appropriate derivatives of (f79*|). that the Hamilton equation 

** - $ («3) 

simply reproduce the results obtained by inverting the definition of the mo- 
menta. 

The action can now be written in Hamiltonian form as 
(p A q A -H(q A ,p A ,N,N i )) d 3 xdt 
= J (P A U ~ NH(q A ,p A ) ~ N t 7i\q A ,p A )) d 3 xdt. (84) 

Variation of this action with respect to p A and q A reproduce the Hamilton 
equations 

«* = w (85) 

and 

P A = -%T> (86) 
dq A 

while variations of shift and lapse yield the superhamiltonian and supermo- 
menta constraint equations 

H(q A ,p A )=0 (87) 
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W(q A ,p A )=0 (88) 

that will be discussed in more detail after we develop convenient notations 
and formal machinery to handle relevant questions. At this point we only 
wish to stress once more that equations (|85|) can be obtained by inverting 
the kinematic relations (J57)l and, thus, can be treated as independent of the 
Hamiltonian dynamics. 



4 Geometro dynamic Superspace. 

York's analysis of the geometrodynamic degrees of freedom suggests that 
the appropriate configuration space for geometrodynamics is not the super- 
space of 3-metrics (or 3-geometries) , but rather the space of conformal 3- 
geometries. We describe here the ideas of such dynamics in a generalized 
form for a case when the 3-metric components are given as functions of 
n q other variables q A , A — 1, . . . , n q < 6, 

9ij = 9iMa)- (89) 

The functions q A are assumed to be independent and form a complete set. 
Following York's analysis, we split the set of variables {q A } n A q =1 into a subset 
{/?r}j=i (nd < 2) of the true dynamic variables and a subset {a M }!|Lo (n e < 3) 
of the so-called embedding variables, with the identifications facilitating a 
comparison with the equations of the previous section: 

Pi = qi (90) 

«/i = Qn d +n+l = Qix- (91) 

This allows us to freely switch notations to better suit the context. It is clear 
that 

n q = n d + n e + 1. (92) 

We wish to reformulate geometrodynamics on the configuration space 
(geometrodynamic superspace) of the true dynamical degrees of freedom qj = 
Pi (conformal superspace as opposed to the superspace of 3-metrics of the 
previous two sections). 

The Lagrangian C, equation (JBJ with K given by (|34j) . remains a function 
of the same variables as in the previous section. It can be written more 
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appropriately as 



£ = C(q A ,q A ,N,Ni) 

= C(qi, qi, g M , g M , N, A 7 *) = £(/3j, /?/, a M , d M , AT, A 7 *), (93) 

where only g/, (g)j are related to the new configuration space. The rest of 
the arguments of £ are functional parameters, which make the Lagrangian 
£ explicitly time dependent (although the time dependence is introduced 
through functions q^, q^, N, Ni). Only the momenta p 1 conjugate to qi 

' = 1 (94) 

retain their dynamical meaning (and eventually become arguments of the 
Hamiltonian). The similar quantities 

*• = £ (95) 

can be introduced and used in developing the theory, but cannot be treated 
as momenta. 

The transition to a Hamiltonian formulation, and the analysis of the 
constraints, can be performed almost the same as in section |21 It is especially 
simple in the case when the matrix Q AB , given by (|56J1 . is such that Q M/ = 
Q 1 ^ = 0. That is, Q AB has the block structure 

(«' 4B ) = ( Q o Q l ) < 96 > 

in which case the inverse matrix Qab has the structure 

(Qab)=( Q " ° ) (97) 



Q 



where Qu is the inverse matrix of Q IJ and Q^ u is the inverse matrix of 
Q^ v . This case is rather common in applications (for instance, it includes all 
diagonal homogeneous cosmologies) . The true meaning is determined by the 
expressions for p 1 and p^, replacing (|57j) with 



/ V9 
p = - — 

y 2N 



-G 13 lm M T ij (N l{m + N l]m ) + Q IJ qj] (98) 
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and 



2N 



G vlrn M H. j{Nllm + Nllm)+ Q»v 4i 



(99) 



It is important to notice that the p^ depend only on time derivatives of q^, 
and do not involve time derivatives of the true dynamical variables. 

All the computations in the previous section can be repeated literally, up 
until the Lagrangian is expressed in terms of momenta, equation (J7SJ). In 
view of ()98|) and ()99|) . the final expression for £ takes the form 



£ 



N 



N 



gi R + g 2 QabP A P b 



gl R + g 2 (Qup 1 p J + Q 



(100) 



which is essentially the same as in the previous section, except now the p M 
are not momenta but merely functions 



p" 



p^QA^^N,^] 



(101) 



given by (J99)) . This makes the Lagrangian £ a function of the new arguments 

£ = C(q A , p 1 , q„, N, = £{q h q, } p 1 , q„ N, N t ). (102) 

The Hamiltonian 1~Ldyn on the dynamical superspace 

Hdyn = p'qi ~ £{q A , P 1 , ? M , N, Ni), (103) 

is distinctly different from the Hamiltonian 7i on the superspace of 3-metrics 
described in the previous section. A useful form of this Hamiltonian can be 
obtained by writing it as 



T-Cdyn = p A qA - £(qA,P*, %, N, -p^q^ 



(104) 



H(q A ,P I ,q l *,N,N i ) 



where the first two terms form TC of the previous section, given by equations 
(JB0J)-(JB2J) , with expressions (|9*9*jl substituted for p^ . More precisely, 



H 



DYN 



Hiq^p^q^N^J-p^ 



(105) 



where 



H{urfAnN,Ni) = NH(qA,p I ,q f i,N,N i )+N l H t (q A ,p I ,q t .,N,N l ) (106) 
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with 



niq^p'A^N.Ni) =g^Q AB p A P B -g*R 

= (Qup i p j + QhvP^p") - a^R 



(107) 



and 



W{qAtf,%,N,Ni) = (-2Q AB p A M B lm G lmi i) ]j 

= -2 (Qjjp 1 M J lm + Q, uP »M» lm ) G^'^108) 

The Hamiltonian action on the geometrodynamic superspace, 



p 1 qi — Hdyn d S x dt 

p'qi - (n(q A , P 1 , N, Ni) - p% 

can be used to derive n d pairs of Hamilton equations 



d xdt 



(109) 



p 



dH 



DYN 



dp 1 

9HdYN 



dqj 



(110) 



by varying p 1 and qi. 

Variations of the action with respect to iV and Ni yield n c constraints, 
with n c < 4. In view of equations (|105|) - (|108|) . the functional dependence 
of Hdyn on N, N can be expressed as 



H DYN {q^P 1 \p"(N, N h qj, N, N 
= NH (q A ,j/,jf{N, N u y) + N t W (q^p^p^N, N h q,) 



(111) 



-jf{N,N h q^q^. 
Variation of N produces the constraint 

OHdyn , OHdyn dpv 



dN 



+ 
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(112) 



H{q A ,p I ,p^N,N i ,q l S), (113) 



and equation (|111|) implies 

ON 

so, together with the trivial non-dynamic expression (|55j) of the previous 
section, 

9TC dyn 97t 



dp- dp- ^ = ff,"^ = 0, (H4) 
which reduces (j!12|) to the superhamiltonian constraint 

H(q A ,p I ,lf(N,N i ,q l S)=0. (115) 

Likewise, variation of iVj leads to the supermomenta constraints 

H t (q A ,p I ,p^N,N l ,q fM ))=0. (116) 

These superhamiltonian and supermomenta constraints on the geometrody- 
namic superspace are obtained from the constraints on the superspace of 
3-geometries by a simple substitution of p- as given by (|99|). 

In the general case the matrix Q AB does not have a block structure. 
Instead, it can be written as 

(o = ( £ ) , di7) 

in which case the inverse matrix Qab has the structure 

It is clear that Qjj ^ Qjj and that Qjj is not the inverse matrix of the 
sub- matrix Q IJ . Likewise, Q^ u ^ Q^ u and Q^ u is not the inverse matrix of 
the sub- matrix Q- v . 

Expressions for p 1 and p M take the form (compare with (J98)) and (|99|l ): 



i \J~9 
p = - — 

y 2N 

and 

vr = — — 
y 2N 



■G^ lm M\ 3 (N l]m + N l]m ) + Q IJ qj + Q Iv q„] (119) 

-G^ lm M\i (N l{m + Ni\ m ) + Q» J qj + q v ] . (120) 
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An important change is that now the p^ depend not only on time derivatives 
of but also on time derivatives of the true dynamical variables. 
Expressing the Lagrangian £ in terms of momenta yields 

C = N [g^R + g-h ( K Q IJ p I p J + Q Iu p T p v + Q,j P »p J + Q^P fl p u )' 

(121) 

which is again essentially the same as in the previous section, except now p^ 
are functions 

p^ — p fJ, {q A - l Qi, N, Ni) (122) 

given by (|12(Jj) . This in turn makes the Lagrangian £ a function of the new 
arguments 

£ = £(q A , p 1 , p»(q A , q r , q„ N, N t ), N, N { ). (123) 

The dependence of the Lagrangian on <jj in the general case seems to break 
the argument developed above in the case when Q AB has block structure. In 
fact, this is not the case, and the problem can be easily remedied. 

As before, we assume that the true dynamical variables qi are indepen- 
dent, which implies that the sub- matrix Q IJ in ()119|) is invertible (we keep the 
notation Qjj for the elements of the matrix inverse to the sub-matrix Q IJ ). 
This means that (|119|) can be considered as a system of linear equations with 
respect to qi. Its solution 

qi = Qu G illm M\ 3 (N l]m - N l]m ) - Qjj Q Jv q v + 2Ng^ Q u p J (124) 

expresses qj as functions 

qi^QiiqAJA^N.Ni). (125) 
Their substitution in (|122j) transforms p^ to 



p" = p» (q A , qi(qA,p T , ? M , N, N t ), q^ N, N>) = p^p 1 , N, N t ) (126) 
and £ into the form 

£ = £ (q^p^p* (u, qiiQAiP 1 , q„, N, N^, q„, N, JV«) , N, N t 

= £(q A ,p I ,q fi ,N,N i ). 



(127) 



With this in mind, we can follow the same steps as before in developing the 
Hamiltonian formalism. We introduce the geometrodynamic Hamiltonian as 

Hdyn = p'qi ~C (128) 



21 



where £ is given by (J121)) . This allows us to recover the same expressions for 
7~Ldyn 

H DYN = p A q A -C- =H- p^, (129) 

with 

H = NH + NiK (130) 

where Tl, 7i, 7i l are given by the same expressions as in the previous section, 
except that p M is now given by ()12()j) and qi is given by ()124j) (also, Qab from 
the previous section should be replaced by Qab)- 

Tracing the chain of arguments developed above results in the action 



d xdt 



I = J ptqj - H DY n 

= J [p'qj- (H-p^q^)] d 3 xdt, (131) 

with TCdyn and TC arguments best described by equations 

Hdyn = T-Ldyn (qA, P 1 , P^ (qa, qi, N, Ni) , N, 

H = U{q A J \f {q A AiA^N,Ni) ,N,Ni) (132) 
qi = qi(qA,P?,qn,N,Ni)- 

Following the same line of arguments as before leads to the Hamilton equa- 
tions ((HE} and to the n c constraint equations 



H^P^P^N^qjiq^p^q^N,^),^)) = (133) 

^{u^^N^NiAMa^A^Ni),^)) = 0. (134) 

These superhamiltonian and supermomenta constraints on the geometrody- 
namic superspace in the general case are obtained from the constraints on 
the superspace of 3-geometries by a simple substitution of p M as given by 
(PUD and qj as given by (IP21) . 



5 The Hamilton— Jacobi Equation. 

A detailed description of the Hamilton- Jacobi equation on the configuration 
space of 3-metrics can be found in the literature, including monographs, 



22 



together with a detailed explanation of how the concept of a functional (bub- 
ble) derivative takes an active part in writing the equation. We therefore 
provide only a very brief discussion of this subject here. 

Following standard approaches we introduce the Hamilton principal func- 
tional, 

S = S[t,g lk ] (135) 

as the extremal value of the action evaluated between the 3-slices given 
by (tf , g r ik ) and (t,gi k )- The primed slice is assumed to be fixed: 

S[t, Qik] = hxtremum = J ^ 9ij ~ ~Hw Dw) d? X dt (136) 

(*'-<4) 

The integral on the right hand side of this expression is extremized with both 
ends fixed. If the upper limit is released, the integral becomes a functional of 
t and given by (|136|) . The Hamilton-Jacobi equation for this functional is 
obtained by variation of the upper limit. Sometimes this variation is preceded 
by imposing the constraints [4|. However, such a move weakens arguments 
following the variational procedure. Instead one can vary the integral before 
imposing the constraints, as it is done in mechanics. It is easy to see that the 
variation (St, Sgik) on the final hypersurface produces variation in S given by 



SS 



J [ir ik Sg lk - H WDW (g ik , vr 4fc ) St] d?x. (137) 



The standard expression for this variation in terms of functional derivatives 
of S is 

.„ f SS SS 
SS = I -j-St+ - — Sg ik 
St Sg ik 



d 3 x, (138) 



and comparison of the two equations yields the expression for momenta 



and one more equation 



n ik = ^- (139) 
Sgik 



5 -^ = -n WDW (g lk ,n ik ). (140) 



Together, these result in a functional differential equation for S 

SS ( 5i 

77 — riWDW Qiki T" 

St \ Sg ik 



SS _ ( 5S\ , 1yl1 , 
— = HwDW 9ik, ~7 — (141) 
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that should be considered the Hamilton-Jacobi equation. However, in Hamil- 
tonian geometrodynamics on the superspace of 3-metrics, the constraints 
and (j3*Uj) change the nature of this equation in such a way that, al- 
though being an important statement concerning the nature of time, it no 
longer performs the same functions as the Hamilton-Jacobi equation in me- 
chanics, at least in the standard treatment of the subject. The standard 
approach is to first impose the constraints, which take the form 

= (142) 

= 0. (143) 

Equation (j2~3j) for H-wdw then allows us to write the Hamilton-Jacobi equa- 
tion (|141|) in the form 

ss . , „ , / ss \ _ _ „ „• / SS \ 

-li =NH {^) +N ' w {^) (144) 

which, together with (fT32jl and (fT33|) yields 

f = 0. (145) 

From the dynamical point of view the partial functional derivative of S with 
respect to t is computed with only g ik fixed. Equation (|145|) might seem 
strange, since (|136|) implies that S also depends on N and iVj, which might 
in turn depend on t. This does not create problems because, as is easy to 
see, 

w-wr - (146) 

This means that equation ()145|) only states that S does not depend explic- 
itly on time and that any dependence on t can emerge only through the 
components of g^. In other words, information about time is carried by the 
3-metric of the slice, which, after proper refinements provides the basis for 
Baierlane, Sharp and Wheeler's concept of intrinsic time. It also implies that 
the functional S[t,gik] is a functional of the slice 3-metric only 

S = S[g ik \. (147) 
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In any case, equation (|14U|) . although reminiscent of the Hamilton-Jacobi 
equation of mechanics, does not perform functions that are expected from 
the Hamilton-Jacobi equation. This statement requires more refined consid- 
erations, based on the observation that equations ()143|) can be written as (cf. 

(f) =0, (148) 

which can be interpreted [3] as invariance of S with respect to the choice 
of coordinates (diffeomorphism invariance), and can be expressed by the 
statement that S is not even the functional of the 3-metric, but only of its 
diffeomorphically invariant part, called the 3-geometry 

S = Sf^]. (149) 

This reduces the left hand side of equation ()144|) . together with the second 
term on the right hand side of the same equation, to the kinematic statement 
expressed by ()149|) . The only remaining part that can possibly have dynamic 
content can be written as 

«( <3 »6.^)-0. (150) 

This equation is identified as the Hamilton-Jacobi equation on the super- 
space of 3-geometries. It can be shown to encode all the dynamic content 
of the theory, but the equation is purely symbolic. In practice, this equation 
can be solved by picking for S[gik] a form satisfying ()148|) . and subsequently 
adjusting the functional parameters of this solution to satisfy (1142}) . Alter- 
natively, we could in principle solve equation (|14Uj) or (|144|) and then adjust 
parameters of the solution to satisfy (|142j) and (J143|) . The first step of this 
procedure produces solutions of ()140|) both on and off the constraint shell, 
which involves providing information (including but not limited to informa- 
tion concerning shift and lapse) that gets eliminated when the solutions are 
restricted to the shell. This procedure is equivalent to the first one, but is 
not as practical in classical geometrodynamics. We will return to it later in 
the context of quantum geometrodynamics. 

The transformation of variables described in sectional and given by equa- 
tions ()31|) . result from replacing the variables go- with the new variables qA- 
This modifies the equations but leaves intact the content of the theory and 
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its interpretation in the generic case. Slight modifications are required in 
non-generic cases, when the dimension of the configuration superspace is 
reduced. 

On the superspace of q A , the Hamilton principal functional 

S = S[t,q A ] (151) 

is 

(t,g A ) 

S[t, q A \ = hxtremum = J (p^QA ~ H(N ', N h q A , p A )) d 3 X dt (152) 

(t',<l' A ) 

where TC is given by equations (J79|) - (|52j). The variational procedure on 
the superspace of q A is similar to the one used on the superspace of g$., and 
yields the expression for the momenta 

P A = ^ (153) 

and the equation 

5 -l = -H(N,N i ,q A ,p A ). (154) 
Together these result in the functional differential equation for S, 

-_ = «^, JV „ 5 ,_j. ( 15 5) 

Following the logic of standard Hamiltonian dynamics this is considered to be 
the Hamilton-Jacobi equation. In view of (|80|h this equation can be written 

as 

-f = ™(- S * («*■©■ < i56 > 

The same line of reasoning as before leads to the constraints 

n(q A ,^-) = (157) 



6S 
8q / 



W (q A , = 0, (158) 
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which can also be expressed as 



m = mr°- (159) 

The three equations (jl56j) - (jl58|) imply, as before, 

f = 0. (1 60) 

The last two equations can be summarized by the statement that S[t, q A ] is, 
in fact, a functional of q A only: 

S = S[q A ). (161) 

Additional considerations, similar to those leading to equations ()148|) - 
(I15()j) . are based on the observation that equations (jlf>8j) can be written as 
(cf. (JH2J)) 

(-2Q AB ^M B lm G lm A =0 (162) 



Sq 



A 



\J 



which can again be interpreted, in the generic case, as invariance of S with 
respect to the choice of coordinates (diffeomorphism invariance). This can 
be expressed by the statement that S is not even a functional of q A , but only 
of the diffeomorphically invariant information carried by q A , the 3-geometry 
(cf. equation (|149|0 . This reduces the left hand side of equation (|156|) . 
together with the second term on the right hand side of the same equation, to 
the kinematic statement expressed by ()149|) . The only remaining part that 
can possibly have dynamic content can be expressed by ()150|) . or in practice, 
by equation (|157j) . which might be identified as the Hamilton-Jacobi equation 
on the configuration superspace of q A . Using (jSTj) . we can also write 

g-iQAB^-¥--giR = o. (163) 

oq A dq B 

All the comments made above concerning the solution of the Hamilton-Jacobi 
equation on the superspace of can be repeated with no essential change 
in the generic case of the superspace of q A . 

The non-generic, degenerate, cases require more attention. A straight- 
forward application of variational principles in these cases might not lead to 
the desired result because of restrictions imposed on variations of the action. 
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These arise from the fixed form of expressions for g^., as well as restrictions 
on variations of the shift and lapse (often referred to as gauge conditions). 

The geometrodynamic superspace has been described in section 0J It can 
be thought of as the configuration space of the true dynamical variables qj. 
The Hamilton principal functional on the geometrodynamic superspace 

S = S[t,qj] (164) 

is given by 

(t,qi) 

S[t,qi] = Iextremum = J (V 4l ~ HdYn) d^xdt (165) 

H DY N = H-p^ (166) 



where 
and as before, 



H = NH + N i H\ (167) 

The arguments of TCdyn and Ti are given by (|132|) . with functions qi and p M 
given by equations ()125|) and ()126|) . In the discussion of the Hamilton-Jacobi 
equation below, it is useful to describe the arguments of Hdy n, T~t, W, and 
W as follows 



H D yn = T-Ldyn [qAiP 1 ,P^ (qA, Qi(qA,p J , %i N > N i)> %i N > N i) > N , Ni 

H = n(q A ,p I , P ^(q A ,q I (q A ,p I ,q^N,N, l ),q^N,N, l ),N,N l ) (168) 

n = H(q A ,p I ,p^(q A ,q I (q A ,p I ,q ll ,N,N i ),q ll ,N,N i ),N } N i ) 

n = W( U ,p I ,pP(uMu,P I ,q»N,Nj,qp,N,N^,N,N^. 

The variational procedure on the geometrodynamic superspace of qi , similar 
to the one used before (variation of the endpoints), yields an expression for 
the momenta conjugate to the true dynamic variables qi, 

V 1 = ^, (169) 
oqi 

and the equation 
S S 

-H D yn {q^p^p" (g A , qiiq^p 1 , q» N, Nj, N, N t ) , N, N, 



St 



(170) 
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which together result in the functional differential equation for S 



SS ( SS ( ( SS \ \ \ 

- — = Hdyn Ua, Ua, Qi Ua, j^, %, N, NA , qtM N, NA , N, nA 

^ ^ (171) 

This is considered to be the Hamilton- Jacobi equation. As before, the con- 
straints are enforced by requirements that can be expressed as the functional 
differential equations 

SS SS 

iN=wr°- (172) 

or equivalently (cf. section 0J 

U (q A , (q A , qj (q A , q„ N, N^j , q„ N, N^j , N, N^j = 0. (173) 

n l (q A , (q A , gi (q A , q„ N, n)j , q„, N, N^j , N, n)j = (174) 

On the geometrodynamic superspace we have, in general, 

f^O (175) 

even on the constraints shell. The reason for this is that, unlike the varia- 
tional derivative of S with respect to t of equation (|16U|) computed at fixed 
q A , the derivative of equation ()175j) is computed at the fixed true dynamic 
variables qi only. 

As before, there are two ways to solve the equations for S, similar to 
those described above for the configuration superspace of q A . The analog of 
the traditional choice would be to solve first the system of equations (|171j) 
and ()174J) . which will partially fix the form of S. After that, solve ([173)1 . 
which is considered the proper Hamilton- Jacobi equation. However, such 
a choice mixes the true dynamical variables with embedding parameters at 
both stages and looses all the advantages of the similar procedure on the 
superspace of all q A . This goes against the structure of the theory and 
implies the loss of any similarity to the Hamilton- Jacobi theory in mechanics. 
The second approach is to solve ()171j) first, considering it the Hamilton- 
Jacobi equation, and then using ()173j) and ()174j) to adjust the functional 
parameters of the solution. The first stage of the procedure provides, in 
principle, solutions of the Hamilton-Jacobi equation both on and off shell, 
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and thus includes information that disappears at the second stage when the 
solution is forced onto the constraint shell. In practice, this procedure is 
more involved than it appears, but in the end is equivalent to the procedure 
on the superspace of qA described above. 

To summarize, the Hamilton-Jacobi theories on the superspace of 3- 
metrics (parametrized by or q^) and on the geometrodynamic superspace 
are found to be equivalent. Any possible difference is erased by forcing solu- 
tions onto the constraint shell. Identification of the Hamilton-Jacobi equa- 
tion is, to an extent, arbitrary. The total content of the Hamilton-Jacobi the- 
ory is expressed by three equations, be it on the superspace of 3-geometries 
(equations (|156j) . (|157|) . and (|158j) ) or on the geometrodynamic superspace 
(equations (|171j) . (|173|) . and (jl74j) V Identification of one of the equations 
as the Hamilton-Jacobi equation is related to the choice of solution strategy 
and does not influence the final solution. It is mostly a matter of convenience 
and, as such, is problem dependent. The situation changes when the theory 
is quantized because it is the Hamilton-Jacobi equation that is converted into 
the wave equation, while the other two equations merely supply additional 
information. 

6 Canonical Quantization 

Canonical quantization of any field theory is based on the Hamilton-Jacobi 
representation of the theory and consists of steps that are determined by this 
representation. Instead of the classical system determined by the functional 
S, a quantum system determined by the state functional \1/ is introduced. 
The arguments of \1/ are assumed to be the same as the arguments of S, 
and are determined by the choice of the configuration superspace and by the 
assignment of the equation to be treated as the Hamilton-Jacobi equation. 

The classical Hamilton-Jacobi equation is then transformed into the wave 
equation of the quantum theory. The functional derivatives of 5* in the 
Hamilton-Jacobi equation are replaced by operators acting on \1/ as follows 



5S 




(176) 




5S 



h 5 



(177) 




h 5 



5S 
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The rest of the expressions participating in the Hamilton- Jacobi equations 
are interpreted as c-numbers acting either on the density ip of ty, defined by 



ty = J ipd d x, (180) 

or the densities produced by functional differentiation of ty. For instance, if 
A is an expression that does not contain momenta, the action of the operator 
A associated to it will be given by 

Aty = Aip (181) 
~ Sty Sty 

A t — = A^- (182) 
Sq A Sq A 

A ir~r- = A ir~r- 183 

Sq A Sq B Sq A Sq B 

and so on. 

On the configuration superspace of 3-metrics, with equation ()142j) identi- 
fied as the Hamilton-Jacobi equation, this procedure produces the Hamilton 
operator (we ignore here technical problems such as factor ordering) 



n = n{ gij ,^) (184) 

and the wave equation 

n( gij ,r j ) ty = o, (185) 

known as the Wheeler-DeWitt equation. 

Similarly, the change of variables that introduces the superspace of q A , 
and the assignment of (|157|) as the Hamilton-Jacobi equation, produces the 
Hamilton operator 

n = n(q A ,p A ) (186) 



and the wave equation 

n(q A ,p A )ty = 0. (187) 



A more detailed form of this (Wheeler-DeWitt) equation, based on (|163|h is 

-n 2 g-^Q AB —- g*R* = 0, (188) 

Sq A Sq B 
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which is a second order functional differential equation. The second func- 
tional derivatives do not participate in the classical Hamilton- Jacobi theory; 
simple descriptions of these operations can be found in the literature jlOj 
(more sophisticated treatments are also easy to find). This equation does 
not contain the time derivative of For it to describe the time evolution of 

time must be inserted into the equation by assigning a suitable function of 
the 3-metric. The resulting equation cannot be interpreted as a Schrodinger 
equation. Its structure is more similar to that of the Klein-Gordon equation, 
and this is the cause of some of the problems of time in the quantum picture 
associated with this approach. 

The wave equation constitutes only a part of the theory. In addition, com- 
mutation relations are imposed on the coordinates of the configuration super- 
space (gij or q A ) and their conjugate momenta (tt^ or p A ). The state func- 
tional is a distribution over the superspace of configurations, which makes 
inserting time into the quantum picture, as described above, more trouble- 
some than it initially appears. 

The quantum version of the auxiliary relations (|143j) and (|158|) (super- 
momentum constraints) are obtained by forming the operators 

7? = H l (g ij , : x ii ) (189) 



or 



n l = n l (q A ,p A ) (190) 

and enforcing the constraints by operator equations that can be written, for 
both cases, as 

7?tf = 0. (191) 

In the case of the configuration superspace of q A this equation is obtained 
from ()162|) . and can be written as (again, ignoring factor ordering issues) 

(-2 - Qab |^ M B lm G lm A =0. (192) 



i Sq A 



\j 



The tendency to enforce the supermomentum constraints as operator equa- 
tions, after the superhamiltonian constraint has been interpreted as the wave 
equation and written as an operator equation, is quite understandable. After 
all, the classical versions of these equations are merely different components 
of an equation that ensures energy-momentum conservation. However, once 
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implemented, these equations generate the problems of time and prevent the 
introduction of any concept of time or of geometrodynamic evolution jH]. 

Comparing this quantum gravity picture with quantum electrodynamics, 
we find that the equations in question play a role similar to that of the 
Lorentz gauge condition (charge conservation). In quantum electrodynamics 
the Lorentz gauge is imposed as a statement concerning expectation values 
rather than as an operator equation[llJ. A similar approach in quantum 
gravity does not seem to be viable in the picture based on the Wheeler- 
DeWitt equation. 

In contrast, canonical quantization on the geometrodynamic superspace 
of qi, with equation (|171|) identified as Hamilton- Jacob i equation, appears to 
produce a quantum gravitational picture that avoids the problems of time, 
and which is similar in spirit to quantum electrodynamics. The geometrody- 
namic quantum Hamiltonian operator is based on the expression for T~Cdyn 
on the right hand side of equation (jl71|) . 

H DY n = T-Ldyn (qA: P^P 1 * (qA: Qi (qa, f ■, %■> N, N^j , g M , N, N^j , N, N-j , 

(193) 

where p 1 is given by 

P 1 = ~ t-. (194) 
i bqi 

The wave equation of the geometrodynamic quantum theory (Schrodinger 
equation), 

Sty ~- 

-ih— = H D Y N y, (195) 
ot 

is obtained from (|17H) in the standard way. 

The commutation relations are imposed only on the true dynamical vari- 
ables qi, and their conjugate momenta p 1 . Embedding variables and their 
velocities are c-numbers, and as such generate only trivial commutation re- 
lations. Time, suitable for describing quantum geometrodynamic evolution, 
can easily be introduced through these embedding variables. 

The auxiliary conditions (|173j) and (|174j) of the classical theory are re- 
placed by quantum conditions, based on operators derived from the left hand 
side of these equations. This is achieved using a procedure similar to the one 
used to form Hdyn' 

H = H{qA,f,P^ (qa, qi (qA,p ! i N, N t ) , q„ N, N t ) , N, N t ) (196) 
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U % = H* (q A , f, {qA, qj (q A , f , q„ N, N { ) , q„ N, N t ) , N, N t ) . (197) 

It has been mentioned above that imposing the constraints in operator 
form leads to numerous difficulties, including time problems. A weaker way 
to impose the constraints, similar to that of quantum electrodynamics, is to 
impose the constraints on expectation values. 

The first step of this procedure is to solve the Schrodinger equation (|195|) . 
with appropriate initial and boundary conditions, assuming that embedding 
variables are represented by c-numbers which are unknown but assumed to 
exist. The resulting solution ty s is a functional that can be represented as 

*s(t,qi) = J Mt,Qi;^)d 3 x. (198) 

The action of any observable A on this produces a function on slices. The 
expectation of this observable over the solution can be written symbolically 

as 

A s = (V s \A\y s ) (199) 
and computed following the prescription 

A s = J r s A^ s V qi Vqj, (200) 

which includes functional integration over the geometrodynamic configura- 
tion superspace but not over slices. 

Quantum constraints on the level of expectations can then be formed as 

(y s \H\^ s ) = (201) 
= 0, (202) 

with H and W given by (115^ and ljI37l) . 

7 Discussion. 

In this paper we have extended the geometrodynamic quantization approach 
from quantum cosmo logical models [6 J to the generic gravitational field. We 
find no inconsistencies in this broader setting. Quantum geometrodynamic 
evolution is determined by the Schrodinger equation ()195|) . together with the 
constraints (|201|) and (|202|) imposed on expectations. Lapse and shift should 
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be either given explicitly or determined by four additional conditions, which 
determine the interpretation of the time parameter t. If the classical version 
of the conditions includes the true dynamical variables and their conjugate 
momenta, their quantum version is imposed in the form of expectations, just 
like the constraints. 

Solving any particular problem can be thought of as a three-step proce- 
dure. First, the Schrodinger equation (|195|) is solved, assuming that the em- 
bedding variables and their time derivatives are unique, although unknown, 
functions of time and the spatial coordinates. The same is assumed of lapse 
and shift, unless they are given explicitly. Solving the the Schrodinger equa- 
tion implies that appropriate boundary and initial conditions for the state 
functional on the geometrodynamic configuration superspace are given. The 
resulting solution is a functional that depends on the embedding variables 
and their time derivatives, as well as shift and lapse. The expectations of 
the constraints (j201|) and ()202j) over the solution of the Schrodinger equation 
are then computed. This procedure produces four differential equations for 
the four embedding variables if the lapse and shift are given explicitly. Al- 
ternatively, one can simply couple the four constraint conditions with four 
functional conditions for the lapse and shift. These procedures determine the 
meaning of time. The last step is to solve these equations and substitute the 
solutions for the embedding variables, their time derivatives, and lapse and 
shift into the expression for the state functional. 

This whole quantization procedure for the general gravitational field par- 
allels the quantum cosmo logical examples considered elsewhere [6, . Consider- 
able complications are caused by the algebraic complexity of the expressions 
for the geometrodynamic Hamiltonian and the constraints, as well as by the 
functional nature of equations. These complications do not, however, stop 
the solution procedure in principle, although they do introduce a rather com- 
plex coupled system. This complexity places demands on our ability to gain a 
proper understanding of the problem, especially in setting appropriate initial 
and boundary conditions on the configuration superspace. 

In principle, however, the geometrodynamic quantization formalism in 
the general setting retains all of the essential features previously illustrated 
in the context of homogeneous cosmologies. 



35 



References 

[1] J. Butterfield and C. J. Isham, "Spacetime and the Philosophical Chal- 
lenge of Quantum Gravity" , |gr"-qc/9903072[ (1999). 

[2] K. Kuchaf, "Canonical Quantum Gravity", in R. J. Gleiser, C. N. 
Kosameh and 0. M. Moreschi, editors, General Relativity and Gravi- 
tation (IOP Publishing, Bristol, 1993). 

[3] C. W. Misner, "Feynman Quantization of General Relativity," 
Rev. Mod. Phys. 29, 497-509 (1957). 

[4] C. W. Misner, K. S. Thorne and J. A. Wheeler, Gravitation (W. H. Free- 
man and Co., San Francisco, 1970). 

[5] Kuchar, K. V.,"Time and Interpretations of Quantum Gravity" in Proc. 
4-th Canadian Conference on General Relativity and Relativistic Astro- 
physics eds. Kunstatter G, Vincent D E, and Williams J G (World Sci- 
entific; Singapore, 1992). 

[6] A. Kheyfets and W. A. Miller, "Quantum Geometrodynamics: 
Quantum- Driven Many-Fingered Time," Phys. Rev. D51, n. 2, 493-501 

[7] A. Kheyfets, D. E. Holz and W. A. Miller, "The Issue of Time Evolution 
in Quantum Gravity," Int. J. Mod. Phys. All, n. 16, 2977-3002, (1996). 

[8] J. W. York, "Role of Conformal Three-Geometry in the Dynamics of 
Gravitation," Phys. Rev. Lett. 28, 1082-1085 (1972). 

[9] M. P. Ryan and L. C. Shepley, Homogeneous Relativistic Cosmologies 
(Princeton Univ. Press, Princeton, NJ, 1975) Ch. 11.4. 

[10] I. M. Gelfand and S. V. Fomin, Calculus of Variations (Dover Pubns, 
Inc., New York, NY, 2000) 

[11] W. Heitler, The Quantum Theory of Radiation (Dover Pubns, Inc., New 
York, NY, 3rd ed., 1984). 



36 



